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ABSTRACT 

(N . 

5— I ' Context. Centaurus A is the closest active galactic nucleus. High resolution imaging using Very Long Baseline Interferometry (VLBI) enables 
us to study the spectral and kinematic behavior of the radio jet-counterjet system on sub-parsec scales, providing essential information for jet 
emission and formation models. 
■ Aims. Our aim is to study the structure and spectral shape of the emission from the central-parsec region of Cen A. 

' Methods. As a target of the Southern Hemisphere VLBI monitoring program TANAMI (Tracking Active Galactic Nuclei with Milliarcsecond 
Interferometry), VLBI observations of Cen A are made regularly at 8.4 and 22.3 GHz with the Australian Long Baseline Array (LBA) and associ- 
ated telescopes in Antarctica, Chile, and South Africa. 

Results. The first dual-frequency images of this source are presented along with the resulting spectral index map. An angular resolution of 0.4 mas 
x 0.7 mas is achieved at 8.4 GHz, corresponding to a linear scale of less than 0.013 pc. Hence, we obtain the highest resolution VLBI image of 
Cen A, comparable to previous space-VLBI observations. By combining with the 22.3 GHz image, which has been taken without contributing 
transoceanic baselines at somewhat lower resolution, we present the corresponding dual-frequency spectral index distribution along the sub-parsec 
scale jet revealing the putative emission regions for recently detected y-rays from the core region by FermifLAT. 
q ' Conclusions. We resolve the innermost structure of the milliarcsecond scale jet and counterjet system of Cen A into discrete components. The 
5_l | simultaneous observations at two frequencies provide the highest resolved spectral index map of an AGN jet allowing us to identify multiple 
, possible sites as the origin of the high energy emission. 
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1. Introduction 

"xT 1 

^ ! At a distance of 3.8 + 0.1 Mpc dHarris et al.ll2010l) . the giant el- 
25 1 liptical galaxy Centaurus A (PKS 1322-428, NGC 5128) is the 
• | closest active galactic nucleus (AGN). Cen A is usually classi- 
""xl" ■ tied a s a Fanaroff-Riley type I radio galaxy (Fanaroff & Rilevl 
' 1 19741) and it exhibits do uble- sided kilopar sec scale jets which 
end in giant radio lobes (IClarke et al.lll992h . These jets are pro- 
. . ■ duced in the vicinity of a s upermassive black hole with a mass 
> ; of M = 5.5 + 3.0 x 1O 7 M (llsraelll 19981: rNeumavetil2010l) . Due 
J^j . to the proximity of Cen A, an angular distance of one milliarc- 
rS 1 second (mas) corresponds to a linear size of just ~0.018pc. This 
\ makes it an exceptionally good laboratory for studying the in- 
nermost regions of AGN, for studying jet formation and colli- 
mation as well as testing different jet emission models. CenA's 
parsec scale structure has been well studied at radio frequencies 
using Southern-Hemisphere Very Long Baseline Interferometry 
(VLBI), revealing a bright jet and a faint counter-jet at a view- 
ing angle of 50-80 degrees to the line of sight, with speeds 
ranging from y = 0-lc to 0.45c measured on d ifferent scales 
dTingavet all 1 19981 120011: lHardcastle et"alll2003l) . On subpar- 
sec scales, the radio jet-counterjet system was re solved with the 
VLBI Space Observatory Program (VSOP) by iHoriuchi etaD 
(120061) . 



Cen A is detected across the spectrum. CGR (9/EGRET found 
the AGN to be a likely source of y-rays (|Hartman et al. 1999) 
and r ecently, the Large Area Telescope (LAT; lAtwood et all 
2009) detector aboard the Fermi Gamma-ray Space Telescope 
(hereafter FermifLAT) detect ed strong y-ray emi ssion from both 
the central (nucle ar) region (lAbdoetalll2010al) as well as the 
giant radio lobes ( Abdo et al. 2010b), while TeV emission was 
detected by H.E.S S. (lAharonian et all2009l) . 

The Tracking Active Galactic Nuclei wi th Austral M illli- 
arcsecond Interferometry program (TANAMI; Oiha et al. 2010£] 
has been monitoring Cen A at two radio frequencies about two 
times a year since 2007 Nov. We present results from the 
TANAMI experiments in 2008 Nov (Sect. [2]!. These included for 
the first time the TIGO and O'Higgins antennas (at 8.4 GHz), as 
well as full participation of the Tidbinbilla 70 m antenna (at 8.4 
and 22.3 GHz) throughout the experiments, resulting in both the 
highest resolution images of Cen A ever made and the highest 
fidelity spectral index map of its mas jet (Sect. [3}. In Sect.|4] we 
discuss the structure of the inner parsec of Cen A and constraints 
on the location of possible y-ray emission regions. 

2. Observations and Data Reduction 

TANAMI observations at 8.4 GHz (2008 Nov 27) and 22.3 GHz 
(2008 Nov 29) were made quasi-simultaneously with the 



* email:cornelia.mueller@ sterawarte.uni-erlangen.de 



1 http : //pulsar . sternwarte . uni-erlangen.de/tanami/ 



2 C. Miiller et al.: Dual-frequency VLBI imaging of Centaurus A at sub-parsec scales 

Table 1. Image parameters 
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8.4 
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2008-11-29 
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1.77 
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2.01 


1.27 


88 


uniform 


8.4 


2008-11-27 


0.45 ± 0.02 
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2.84 


3.2 ±0.1 


0.68 


0.43 


33 
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2008-11-29 
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1.46 


0.19 


3.3 ±0.1 


1.55 


1.21 


-75 



Columns show: (1) weighting scheme; (2) observing frequency; (3) observation epoch; (4) RMS noise in the clean images; (5) peak flux density; 
(6) brightness temperature derived from T B = 1.22 • 10 12 (S„/[Jy])(v/[Hz])~ 2 (#/[mas])~ 2 using the restoring beam dCondon & Mitchell 1982J); (7) 
total integrated flux density; (8,9,10) major axis, minor axis, and position angle of the restoring beam 
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Fig. 1. O,v)-coverage at 8.4 GHz (grey) and at 22.3 GHz (black) 
for CenA (see Table [TJ. Note that the transoceanic antennas 
O'Higgins and TIGO do not support observations at 22.3 GHz. 

Australian Long Baseline Array (LBA) and associated tele- 
scopes, the 9 m German Antarctic Receiving Station (GARS) 
in O'Higgins, Antarctica, the 6 m Transportable Integrated 
Geodetic Observatory (TIGO) in Chile, both of which are oper- 
ated by the Bundesamt fur Kartographie und Geodasie (BKG), 
Germany, and the 70 m DSS43 antenna of the NASA Deep 
Space Network in Tidbinbilla, Australia. The LBA consists of 
five telescopes, Parkes (64 m), ATCA (5x22 m), Mopra (22 m) 
all in New South Wales, Hobart (26 m) in Tasmania, and Ceduna 
(30m) in South Australia. Figure [TJ shows the resulting (u,v)- 
coverage for Cen A at both frequencies. Note that the resolution 
is lower at the higher frequency as GARS and TIGO cannot ob- 
serve at 22.3 GHz. Data were recorded on LBA Dis k Recorders 
and c orrelated on the DiFX software correlator dDeller et al.l 
120071) at Curtin University in Pert h, Western Austra lia and cali- 
brated and imaged as described in lOiha et all (1201 Ol) . 

3. Results 

Figure [2] shows images of CenA at 8.4 GHz and 22.3 GHz. 
Image parameters and observation characteristics are listed in 
Table [TJ Using uniform weighting at 8.4 GHz, an angular res- 
olution of 0.7 x 0.4 mas at a position angle (P.A.) of ~33°) is 
achieved corresponding to linear scales of less than 3500 AU. 

The most important constraint on the fidelity of the 8.4 GHz 
image is the lack of intermediate baselines between the LBA 
and the two transoceanic ant ennas, evident fro m "holes" in the 
(w,v)-plane (Fig. [JJ; see, e.g.. lLister et al.ll2000i for a discussion 
of possible effects due to uneven coverage of the («,v)-plane). 
To test the reliability of the imaging, we continuously checked 
whether our clean model describes the uncalibrated visibility 



data consistently. In addition, we gauged the robustness of fea- 
tures in given regions of the map by iteratively testing the effect 
of self-calibrating the data without including these features in 
our clean models. We estimate the calibration uncertainties to 
be on the order of 15 % at both frequencies. On-source errors 
due to the uneven (M, y)-coverage were measur ed following the 
methods described bv lTingav & Murphvl(l2001l) and found to be 
smaller than 15% along the inner ~ 15 mas of the jet. Beyond 
this, the on-source errors gradually grow and can exceed 100 % 
towards the counterjet at 8.4 GHz. We also compared our final 
images with previous and su bsequent, somewha t lower resolu- 
tion, TANAMI observations (Miil ler et"aI1 1201 Oh by identifying 
the same main bright features to achieve a self-consistent multi- 
epoch kinematics model which will be analyzed in a subsequent 
paper. 

At both frequencies, our images show a well collimated jet 
with an opening angle 6 < 12° at a mean P.A. of 9 ~5Q° and 
a fainter counterjet (9 — 130°) with an emission gap in be- 
tween. At 22.3 GHz, the counterjet is considerably weaker: for 
the 8.4 GHz counterjet feature around position R.A. = -18, 
Dec = -14, we put a limit of the spectral index a < -2.5 
(defined via F v oc v +a ). Both the P.A. and the opening angle are 
remarkably similar to the corresponding values of P.A. arcm i n = 
50°-60° and flarcmin ~1 2° seen on scales o f < 4arcmin (4kpc) 
from the central engine (Clar ke et al.l!992b . The feature with the 
highest peak-flux density (see Table [TJ has an inverted spectrum 
(see below) and therefore shares the characteristic properties of 
the compact feature at the upstream-end of most blazar jets, typ- 
ically called "the core" of the jet. Note that in the case of Cen A 
our linear resolution is orders of magnitude better than in most 
blazar-jet images and consequently most of the structure seen 
here within the central parsec would be contained in "the core" 
if Cen A was observed at a smaller angle to the line of sight or at 
larger distance. The significant emission features within the sub- 
parsec scale jet seen in the 2007 November 8.4 GHz TANAMI 
image (Miille r etal.ll2010t lOiha et alJl2010h are in good agree- 
ment with those in this image observed one year later. A pos- 
sible widening and subsequent narrowing of the jet appears at 
~25 mas (« 0.45 pc), and also in the earlier lower resolution im- 
age, as well as in the 22.3 GHz image (see Fig.[2]i. The peak-flux 
densities at 8.4 GHz are comparable, 0.6 Jy in 2007 November 
and 0.7 Jy in 2008 November, indicating only a moderate radio 
flux variability. The 2008 November image reveals a small coun- 
terjet displacement from the jet axis. The overall jet structure is 
also s imilar to the earlier 5 GHz VSOP image (Horiuc hi et al.l 
2006). This similarity suggests that at least some characteristic 
structures are moving slowly compar ed to the jet speed mea- 
sured in earlier work s (0.1c-0.45c: iTingav et all I1998L I200TL 
lHardcastle et"aTll2003l) . 

The spectral distribution along the jet constrains the physi- 
cal properties of the radio emission. The synthesis imaging self- 
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Fig. 2. 22.3 GHz and 8.4 GHz images of Cen A from 2008 November 27 and 29, using natural weighting. The inset shows the inner 
region of Cen A at 8.4GHz using uniform weighting. The lowest contours denote the 3<x-noise-level (see TableQ]). All axes are in 
units of mas relative to the phase center. The ellipse in the lower left corner of each panel indicates the restoring beam. 



calibration procedure leads to the loss of the absolute source po- 
sition. Therefore, the frequency-dependent VL BI shift of self - 
absorbed, optically thick emission features (e.g.. lLobanovll 998) 
has to be taken into account to assemble a spectral index map. 
We determined the offset by aligning the positions of optically 
thin Gaussian model components, which are not affected by the 
shift, to be Aa K \ = -0.25 mas and A6 K \ = -0.2 mas. 



4. Discussion and Conclusions 

Recently, Fermif LAT has detected y -ray emission from the core 
region of Cen A dAbdo et al.ll2010al) but the LAT positional ac- 
curacy does not allow to locate the y-ray emission to an ac- 
curacy of better than ~0°1. Our simultaneous dual-frequency 
VLBI data allow us to study spectral changes at sub-parsec 
scale and therefore to constrain the possible origin of the y- 
ray emission from the core region of Cen A. To do this, we fit 
Gaussian model components to the visibility data (see Table |2). 
We can clearly resol ve the 7 mas large core region imaged by 
iHoriuchi et al.1 (120061) into two separate jet features. The brighter 
feature at the base of the jet, with an angular extent of ~3 mas, 



can be described with three Gaussian model components at the 
higher-resolution 8.4 GHz image while two components are suf- 
ficient to represent the lower-resolution 22.3 GHz core region. 
We identify this region as the innermost visible part of the jet 
at this frequency (more central parts being self-absorbed). The 
second bright emission feature at ~3.5 mas downstream at both 
frequencies was also detected by TANAMI in 2007 November 
(Oiha et al.l2010h and might corres pond to the stationary compo- 
nent detected by Tingav et al.l d200ll) at ~4 mas. Co-spatial with 
the jet- widening region at -23 mas at 8.4 GHz, the 22.3 GHz 
image suggests a possible emission component, which might be 
interpreted as an absorbing feature. However, the limited (u, v)- 
coverage does not allow us to fully constrain the exact spatial 
distribution of this emission feature and consequently spatially 
resolved spectroscopic analysis has to be treated with caution. 

Figure |3]represents the first sub-parsec scale dual-frequency 
image of Cen A. The core region has an inverted spectrum (spec- 
tral index < a < 2) which changes from flat to steep down- 



stream. Remarkably, the spectrum remains flat for the inner few 
milliarcseconds of the jet, indicating multiple optically thick 
emission regions in the core region although a gradual steep- 
ening behaviour can be seen. These optically thick regions can 
be identified as multiple possible production sites of the high 
energy photons because they all show high brightness temper- 
atures, compact, structures, and flat radio spectra. This should 
be considered in the discussion of single-zone versus multi- 
zone emission models of the broadband spectral energy distri- 
bution of Cen A and other extragalactic jets. Recent stud ies of 
Cen A (lAbdo et al.ll2010atlRoustazadeh & Bottcherl l20TTl) mod- 
eling the full radio to TeV range of the broadband spectrum have 
used TANAMI spectral data from the inner 0.9 mas x 0.29 mafl 

As Fig.[3]shows, the peak of the innermost counterjet feature 
at 22.3 GHz is offset by about 1 mas to the south-west with re- 
spect to the 8 GHz image. This offset could be interpreted as a 
sign of free-free absorption in an obscurin g torus as observe d 
in other two-sided extragalactic jets (e.g., Kadler et"ail 12004). 
In this case, however, it is puzzling that the well-defined 8 GHz 
component at Aa = -4 mas, A6 = -2 mas seems to have no ob- 
vious counterpart at 22.3 GHz. Free-free abso rption has been ob- 
serve d before in the central region of Cen A dTingav & Murphv 
2001 ). In fact, this may also affect the innermost part of the east- 
ern jet and could be an alternative explanation for the observed 
flat spectral index distri bution. 

In contrast to M87 dJunor et al.ll 19991) . no rapid broadening 
of the jet opening angle is seen, although the linear resolutions 
are comparable. The jet of Cen A is therefore already well col- 
limated on scales of 0.1 pc = 10 15 m. Beyond this distance, the 
opening an gle is ~12°, abou t the same value as measured on 
kpc-scales dClarke et al.ll 19921) . However, the mass of the central 
black hole of M87 is two orders of magnitude higher than that 
of Cen A so we are accessing a very different part of the jet as 
measured in Schwarzschild radii. 

We plan to determine the evolution of spectral index and 
the proper motion using observations from multiple TANAMI 



2 Note the slightly different model com ponent values as compared to 
the values reported in Abdo et al. (2010a), which reflect self-calibration 
and model-fitting uncertainties. 
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Fig. 3. Top: Spectral index map derived from the flux densities at 
8.4 GHz and 22.3 GHz. The flux density distribution at 8.4 GHz 
is shown in solid black contours. The dashed blue contour shows 
the 3<x-noise level at 22.3 GHz. Both maps have been restored 
with a beam of 1.61 x 1.02 mas (6 = 88°) represented by the 
gray ellipse in the lower left corner. Middle: Flux density distri- 
bution at 8 GHz (red) and 22.3 GHz (blue) for a narrow strip 
of 0.1 mas along the jet at a P. A. =50° corrected for the shift 
of Aa = -0.25 mas and A6 = -0.2 mas. Bottom: Spectral in- 
dex given by the flux densities at P. A. =50°. Uncertainties corre- 
sponding to the absolute calibration uncertainties and on-source 
errors of ~15 % and considering the image rms at both frequen- 
cies are shown as shaded regions around the best-fit distributions 
(solid lines). 

epochs. Combining this with y-ray monitoring, this approach 
will allow us to further constrain the possible y-ray emission re- 
gions. Future ultra-high resolution multi-frequency observations 
with a global VLBI array (including the southernmost VLBA an- 
tennas) will let us measure the spectral turnover frequency and 
probe the presence of a dense free-free absorber. In combination 
with kinematic data from the TANAMI observations, we will be 
able to estimate the magnetic field and relativistic electron en- 
ergy density, key parameters of jet broadband emission models. 



Table 2. Modelfit parameters for core region (natural weighting) 
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